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We investigate the quantum phase transitions of the extended Hubbard model at half-filling 
with periodic boundary conditions employing the entanglement of particles, as opposed to the more 
traditional entanglement of modes. Our results show that the entanglement has either discontinuities 
or local minima at the critical points. We associate the discontinuities to first order transitions, and 
the minima to second order ones. Thus we show that the entanglement of particles can be used to 
derive the phase diagram, except for the subtle transitions between the phases SDW-BOW, and the 
superconductor phases TS-SS. 
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I. INTRODUCTION 

The connection between two important disciplines of 
Physics, namely quantum information theory and con¬ 
densed matter physics, has been the subject of great in¬ 
terest recently, generating much activity at the border 
of these fields, with numerous interesting questions ad¬ 
dressed so far [T|. In particular, the properties of entan¬ 
glement in many-body systems, and the analysis of its 
behavior in critical systems deserve special attention. 

In this work we deal with the entanglement of indis¬ 
tinguishable fermionic particles in the one dimensional 
extended Hubbard model (EHM). We focus in the half- 
filling case. The model is a generalisation of the Hubbard 
model El El, which encompasses more general interac¬ 
tions between the fermionic particles, such as an inter¬ 
site interaction, thus describing more general phenomena 
and a richer phase diagram. Precisely, it is given by, 
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where L is the lattice size, aj a and a^ a are creation 
and annihilation operators, respectively, of a fermion 
with spin a at site j, h^ G = oj a ctj,ai 
and we consider periodic boundary conditions (PBC), 
L + 1 = 1. The hopping (tunnelling) between neighbor 
sites is parametrized by t , while the on-site and inter¬ 
site interactions are given by U and V, respectively. De¬ 
spite the apparent simplicity of the model, it exhibits 
a very rich phase diagram, which includes several dis¬ 
tinct phases, namely: charge-density wave (CDW), spin- 
density wave (SDW), phase separation (PS), singlet (SS) 
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and triplet (TS) superconductors, and a controversial 
bond-order wave (BOW). A more detailed description of 
the model and its phases will be given in the next section. 

Our numerical analysis is performed employing entan¬ 
glement measures for indistinguishable particles intro¬ 
duced recently HH, in conjunction with the density- 
matrix renormalisation group approach (DRMG)f7j [8], 
which has established itself as a leading method for the 
simulation of one dimensional strongly correlated quan¬ 
tum lattice systems. DMRG is a numerical algorithm for 
the efficient truncation of the Hilbert space of strongly 
correlated quantum systems based on a rather general 
decimation prescription. The algorithm has achieved un¬ 
precedented precision in the description of static, dy¬ 
namic and thermodynamic properties of one dimensional 
quantum systems, quickly becoming the method of choice 
for numerical studies. 

The paper is organised as follows. In Sec. |TT] we re¬ 
view the model and its phase diagram. In Sec. |III| we 
present the distinct definitions of entanglement in sys¬ 
tems of indistinguishable particles, focusing on the notion 
of “entanglement of particles”. In Sec. |IV| we present our 
results. We conclude in Sec. m 


II. EXTENDED HUBBARD MODEL 

In this section we give a detailed description of the 
extended Hubbard model P0, and its distinct phases in 
the half-filling case. The reader familiar with the subject 
may skip this section. 

Many efforts have been devoted to the investigation of 
the EHM’s phase diagram at half-filling, using both ana¬ 
lytical and numerical methods png. Despite the appar¬ 
ent simplicity of the model, it exhibits a very rich phase 
diagram which includes several distinct phases: charge- 
density wave (CDW), spin-density wave (SDW), phase 
separation (PS), singlet (SS) and triplet (TS) supercon¬ 
ductors, and a controversial bond-order wave (BOW). 
See Fig for a schematic drawing of the phase diagram 
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FIG. 1: (Color online) Phase diagram of the half-filled ex¬ 
tended Hubbard model in one dimension. The distinct phases 
correspond to the charge-density wave (CDW), the spin- 
density wave (SDW), phase separation (PS), singlet (SS) and 
triplet (TS) superconducting phases, and bond-order wave 
(BOW). The order of the quantum phase transitions is iden¬ 
tified by the different line shapes. The order of the two su¬ 
perconducting phases transition (blue dotted line) is contro¬ 
versial, being identified as a BKT transition m, or a second 
order transition f§]. 


at half-filling. 

In the strong coupling limit (|C/|, | V"| t ), one can 

qualitatively characterize its phases as given by a charge- 
density wave, spin density wave and a phase separa¬ 
tion. For a strong repulsive on-site interaction (U > 0, 
U V), the ground state avoids double occupancy and 
the spin density is periodic along the lattice, leading to 
an antiferromagnetic ordering, namely spin-density wave. 
Its order parameter is given by, 

O sdw (k) = \ £ e ifc(TO -" } [«<> - «> «>], (2) 
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where cr| = } ft — t) ■ I n the limit U —> oo, the ground 
state is dominated by the following configurations: 
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where the state is described in the real space mode rep¬ 
resentation, in which each site can be in the following 
set of configurations: {|0), |t), |i), |tl)}- Considering a 
strong repulsive inter-site interaction [V > 0, V U), 
a periodic fermionic density is generated, leading to a 
charge-density wave. Its order parameter is given by, 

O cdw (k) = l Y, l^rnflu) - (fl m ) (flu)} ■ (4) 
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In the limit V —> oo, the ground state is dominated by 


the following configurations, 

(in,o,n,o,...,nL-i,o> + 

+ |o,n,o,n,-,o,UL>). (5) 


In the range of strong attractive interactions (U, V < 0 
or U > 0, V < 0 with |V| S> |f7|), the fermions cluster 
together, and the ground state becomes inhomogeneous, 
with different average charge densities in its distinct spa¬ 
tial regions. Such a phase is called phase separated state. 
In the limit V —> —oo, the ground state is dominated by 
the following configurations, 
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n,n,-,n ( |),o,...,o), (6) 


where {11} is the set of translation operators. 

In the weak coupling limit, different phases appear. 
For small attractive inter-site interactions (V < 0), su¬ 
perconducting phases are raised, characterized by the 
pairing correlations, 


A r = 


E a i>T°j+*4’ 
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with the respective order parameter O s = 
Ex,*' ( A £ A *')- If the on-site interactions are lower 
than the inter-site interactions (U < 2V), the fermions 
will pair as a singlet superconductor, characterized 
by nearest-neighbor (A SSrin ) or on-site (A SSo ) singlet 
pairing correlations given by, 
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where x = 1. On the other hand, if the on-site inter¬ 
actions are higher than the inter-site interactions (U > 
2V), we have a triplet superconductor, characterized 
by nearest-neighbor triplet pairing correlations (A tSriii ) 
given by, 
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where x = 1. 

Note that the difference between the singlet and triplet 
pairing correlations is simply a plus or minus sign. It 
can be clarified if we consider, for example, the case of 
two fermions in a singlet or triplet spin state, given by 
(| ij) ± | ji)) (|ti) =F Expanding this state, we have, 


I *j) (Iti) =F lit)) ± I ji) (Iti) =F lit)) 

= |* to i) T I* i.J t) ± I j t,U) - I j i,t t) 
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( 11 ) 
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where the lower/upper sign corresponds to the sin¬ 
glet/triplet pairing correlation. 

The last phase in the diagram is the controversial 
bond-order-wave (BOW). By studying the EHM ground 
state broken symmetries, using level crossings in excita¬ 
tion spectra, obtained by exact diagonalization, Naka¬ 
mura m argued for the existence of a novel bond-order- 
wave phase for small to intermediate values of positive U 
and V, in a narrow strip between CDW and SDW phases. 
This phase exhibits alternating strengths of the expecta¬ 
tion value of the kinetic energy operator on the bonds, 
and is characterized by the following order parameter, 

O bow (k) = I [{B m , m+1 B n , n+1 ) 

m,n 

— (Bm,m+ 1) (Bfi,n+l)] > (12) 

where = J2a( a ln,er a m+i,a + H.c.) is the kinetic 

energy operator associated with the mth bond. Naka¬ 
mura argued that the CDW-SDW transition is replaced 
by two separate transitions, namely: (i) a continuous 
CDW-BOW transition; and (ii) a Berezinskii-Kosterlitz- 
Thouless (BKT) spin-gap transition from BOW to SDW. 
Such remarkable proposal was later confirmed by several 
works EM], employing different numerical methods, 
like DMRG, Monte Carlo or exact diagonalization. Nev¬ 
ertheless, while the BOW-CDW phase boundary can be 
well resolved, since it involves a standard second order 
(continuous) phase transition, the SDW-BOW boundary 
is more difficult to locate, for it involves a BKT transi¬ 
tion in which the spin gap opens exponentially slowly as 
one enters the BOW phase. The precise location of the 
BOW phase is then still a subject of debate. To the best 
of our knowledge, the best estimates for the transitions, 
taking U/t = 4, correspond to a CDW-BOW transition 
at V/t « 2.16 [TT][T2JE][T2], and to a BOW-SDW tran¬ 
sition in the range V/t « 1.88 — 2.00 [111 IT21 [T51 HR] , or 
V/t = 2.08 ± 0.02 El- 


III. ENTANGLEMENT OF 
INDISTINGUISHABLE PARTICLES 

Despite widely studied in systems of distinguishable 
particles, entanglement or more general notions of quan¬ 
tum correlations have received less attention in the case 
of indistinguishable particles. In this case, the space of 
quantum states is restricted to symmetric (5) or anti¬ 
symmetric (A) subspaces, depending on the bosonic or 
fermionic nature of the system, and the particles are no 
longer accessible individually, thus eliminating the usual 
notions of separability and local measurements, and mak¬ 
ing the analysis of correlations much subtler. In fact, 
there are a multitude of distinct approaches and an on¬ 
going debate around the entanglement in these systems 
[6J HMD]. Nevertheless, despite the variety, the ap¬ 
proaches consist essentially in the analysis of correlations 
under two different aspects: the correlations genuinely 


arising from the entanglement between the particles (en¬ 
tanglement of particles) [0 ITO^f] , and the correlations 
arising from the entanglement between the modes of the 
system (entanglement of modes) |25lB5] . These two no¬ 
tions of entanglement are complementary, and the use of 
one or the other depends on the particular situation un¬ 
der scrutiny. For example, the correlations in eigenstates 
of a many-body Hamiltonian could be more naturally 
described by entanglement of particles, whereas certain 
quantum information protocols could prompt a descrip¬ 
tion in terms of entanglement of modes. Once one has 
opted for a certain notion of entanglement, there are in¬ 
teresting methods to quantify it |[H0 HTTl HTT] . 

Entanglement of modes can be understood by mapping 
the quantum state in its number representation, namely, 

4-4 \ vac ) —> io...i n ...i 7JV ...o) 

(13) 

where j- L = 1,..., M, and { a]}jLi is an arbitrary set of M 
fermionic operators describing the single particle modes 
of the system (not necessarily the real space modes as in 
the Hamiltonian definition). We will denote hereafter as 
“configuration representation (number representation)” 
the left (right) side of the previous equation. Such equa¬ 
tion corresponds to a mapping to distinguishable qubits, 
represented by the occupied (|1) • ) or unoccupied (|0) ? ) 
modes, which then allows one to employ all the tools 
commonly used in distinguishable quantum systems in 
order to analyze their correlations. One could, for exam¬ 
ple, use the von Neumann entropy of the reduced state 
representing a block with i modes, in order to quantify 
the entanglement between this block with the rest of the 
modes. The reduced state is obtained by the partial trace 
in the number representation (pe = Tr^i (\ip)(ip\))■ No¬ 
tice that, in the mode representation, local observables 
may actually involve correlations between particles. For 
example, in the Hubbard model, although the operator 

f.” acts locally at the jth site (real space 
modes), it describes pairing correlations between par¬ 
ticles. The algebra of local observables at the modes, 
defined in the number representation, is generated by, 

D; oc = |(9i ® J 2 ,m ; 2 g> 0 2 ® 2 3 , m ;■■■ 

••• ;2 1j(m _i)®O m II 4 = 4} (14) 

where 2=2 \ ® Ii+\ (g) • • • (g> 2with j > i and 2 \ 
is the identity operator acting on mode i. In this way, 
unentangled states are those which can be completely 
described by such local observables. It is known that 
such states are simply the separable states in the usual 
tensor product form, \ip) un = \4>i) <g> ^ 2 ) •••<£> \4>m)- 

Based on the previous reasoning, we now define the 
notion of entanglement of particles. Notice first that 
one cannot analyze the system under the usual paradigm 
of separability and locality , where the reduced states ob¬ 
tained by partial trace are mixed ( p r = Tr 2 . ..Jv(lV’XV’l)), 
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whenever the global state is pure and entangled. There¬ 
fore, in the case of indistinguishable particles in the con¬ 
figuration representation, the use of partial trace to char¬ 
acterize entanglement should be carefully reviewed, since 
it would suggest that all pure fermionic states are entan¬ 
gled, given that their reduced states are always mixed. 
In order to generalize the notion of entanglement for sys¬ 
tems of indistinguishable particles, the approach based 
on the algebra of observables sheds light on the problem 
and allows us to go beyond the paradigm of separability 
and locality. 

We now define the proper algebra of “local observables” 
as the one composed by operators which do not cre¬ 
ate correlations between the indistinguishable particles. 
Such algebra, defined in the configuration representation, 
is generated by the following single particle operators, 

Qioc = |o <B> 1 2 ,n + X ® O ® X 3 ,jv + • • • 

• • • + l) ® O || = O j , (15) 

where N is the number of particles. Equivalently, using 
the second quantization formalism, the above set is given 
by the number conserving quadratic operators, f2; oc = 
{(a\aj + H.c.) | i,j = 1, ...,M}. The states that can be 
completely described by such algebra form, in this way, 
the set of unentangled states, where any particle is not 
entangled with any other. Intuitively, we would expect 
that this set corresponded to single Slater determinants 
with fixed particle number. More precisely, for a system 
with N fermions, it is given by, 

Wun = a h a ] 2 - a l N \ VaC ) > ( 16 ) 

where {a]} is an arbitrary set of fermionic operators. Re¬ 
call that these operators cannot be quasiparticles with 
particle-hole superpositions, as usual in a Bogoliubov 
transformation, since the above states have a fixed num¬ 
ber of fermions. In fact, distinct approaches confirmed 
that such set does indeed correspond to the unentan¬ 
gled states E2 QZH23I- The only non-classical correla¬ 
tion present in such states is the exchange , due to the 
antissymetrization, which does not constitute entangle¬ 
ment. For example, in m the analysis follows by us¬ 
ing a very elegant mathematical formalism, called GNS 
(Gelfand-Naimark-Segal) construction, for the case of 
two fermions, each one with Hilbert space dimension 3 or 
4, and two bosons with dimension 3; in ESI Ell the au¬ 
thors propose a “Generalized Entanglement (GE)” mea¬ 
sure, obtaining a simple formula for the “partial trace”, 
and the set of fermionic unentangled states for an arbi¬ 
trary number of particles; or also in f6], where a general 
notion of quantum correlation beyond entanglement (the 
quantumness of correlations) is investigated by means of 
an “activation protocol”, which yields the same set of 
states with no quantumness as the above unentangled 
one. 

As in the case of distinguishable modes, the von Neu¬ 
mann entropy also provides a good quantifier for the en¬ 


tanglement of indistinguishable particles. We can define 
the Shifted von Neumann entropy of entanglement (4j as 
follows, 


E p (\il>M) = S(p r )-]og 2 N, (17) 

where p r = Tri...TrN-i(\ip)(ip\) is the single particle re¬ 
duced state, the partial trace is taken in the configura¬ 
tion space, and S(p ) = Tr(—p log 2 p) is the von Neumann 
entropy. Such a quantifier is obtained simply by notic¬ 
ing that each extremal state in the one particle reduced 
space is respective to a unique single Slater determinant 
[55] . More precisely, the single particle reduced state of 
a single Slater determinant (as Eq. ©) is given by, 

1 N 

^ a). | vac) (vac\ a j( (18) 

V i —1 

and its particle entanglement is null, EpQtj)) (ip\) = 0. 
If the state cannot be described by a single Slater 
determinant, its entanglement is necessarily non null, 
E p ( ||) > 0, and at least one of its particles is en¬ 
tangled with another one. Maximally entangled states 
have their single particle reduced states described by the 
maximally mixed state, as for example, the strong cou¬ 
pling limit phases SDW, CDW, and PS, as described 
in Eqs. (§, (§, (§, respectively, whose reduced state is 
given by, 

1 L 

Pr= 7^- 5Z \ VaC ) ( vac \ ( 19 ) 

i=l <7=1-4 


which have maximal von Neaumann entropy, S(p r ) = 
logn(2 L), and consequently maximal particle entangle¬ 
ment, E p (\ip){ip\) = 1. 

If the Hamiltonian has certain symmetries, its ground 
state entanglement can be analytically calculated as a 
simple function of its quadratures @j. In our particular 
case, from both S z and translational symmetries in the 
extended Hubbard model, formally given by, 

p9 ) = °’ ( 2 °) 

<77^ a 

Tr ( a \a a ja P 9 ) = Tr(a\ i+SyT a {j+S)a p g ), (21) 


where is the annihilation (creation) fermionic op¬ 
erator of a particle in the jth site, with spin ex, and 
Pg = is the ground state of the Hamiltonian, we 

have that its single particle reduced state (p r (ia, ja) = 
ji Tr(aj & ai a pg)) is disjoint in the subspaces with dis¬ 
tinct spin, p r = 0 p^=4 ) an d each of these terms is 

given by a circulant matrix, 


/ X 0 X 1 • • • Xl-2 X L -i\ 

XL—! Xo x 1 Xl-2 
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x 2 

\ Xi 


XL -1 Xo 


X 2 


X\ 
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(23) 


where L is the lattice size. In this way the matrix is 
easily diagonalized, and its eigenvalues {A£} are given 
by a Fourier transform of the quadratures, 


L—l 


XI = k = 


6 =0 


„ 27T 27r 

0, —,(i 1) — 


(24) 


The entanglement is then directly obtained from Eq. (SH- 


iv. ENTANGLEMENT AND QUANTUM 
PHASE TRANSITIONS 

The computation of the single particle correlations, 
and consequently the entanglement of particles, was nu¬ 
merically performed using DMRG. Although DMRG is 
less accurate for problems with periodic boundary condi¬ 
tions (PBC) than with open boundary conditions (OBC), 
from the Physical viewpoint PBC are strongly preferable 
over OBC, as boundary effects are eliminated and finite 
size extrapolations can be performed for much smaller 
system sizes. In this work we analyze the extended Hub¬ 
bard model considering PBC. Our simulations were per¬ 
formed for systems up to L = 352 sites, always keeping 
a large enough dimension (m) for the renormalized ma¬ 
trices (ranging from m = 100 to 1000) and number of 
sweeps (~ 20 sweeps), in order to obtain an accurate 
precision. In Fig(2] we see that m ranging from 200 to 
300 is enough for an entanglement accuracy of the order 
of 0( 10” 4 ). The accuracy for the ground state energy, as 
well as the truncation error, using such parameters, are 
of the order of G(10~') 

Our results for the entanglement of particles in the ex¬ 
tended Hubbard model at half-filling are shown in Figj3] 
to be dissected below. It is remarkable that such picture 
highlights the known phase diagram of the model. We 
first note that, as expected, we have a maximum of entan¬ 
glement at the strong coupling limits (E p —> 1), and as we 
decrease the interactions between the particles, the en¬ 
tanglement tends also to decrease, until the unentangled 
case for the non interacting Hamiltonian (U = V = 0). 
The figure thus presents the shape of a valley around 
this point. Following then the discontinuities and the lo¬ 
cal minimum points in the entanglement, we can easily 
identify the quantum phase transitions, except for both 
the subtle SDW-BOW transition, and the transition be¬ 
tween the superconductor phases TS-SS. In the former 
case, one needs to recall that the observation of the BOW 
phase is by itself a hard task, since its gap opens expo¬ 
nentially slowly, and also that there are evidences that 
such transition is of infinite order 123133 ■ Therefore we 
believe that a possible detection of such transition by the 
entanglement of particles would require higher precision 
numerical analysis as well as the study of larger lattice 
sizes. Concerning the TS-SS transition, on the one hand 
the order of the two superconducting phases transition 
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FIG. 2: (Color online) Accuracy analysis for the computa¬ 
tion of entanglement of particles using DMRG. It is shown the 
accuracy of the entanglement, A (E p ) = E p (m) — E p (m — 50), 
as a function of m (dimension of the renormalized matrices), 
at the point U = 4, V = 2.11, and using 20 sweeps in the com¬ 
putation, which is enough for the ground state convergence. 
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is controversial, being identified as a BKT transition m 
as well as a second order one [H] in the literature. On 
the one hand, we would be led to strengthen the result 
of a BKT transition, since our entanglement does not 
detect it. On the other hand, it is reasonable the apa¬ 
thy of the entanglement of particles on distinguishing the 
two phases, since the correlations between the particles 
in the two superconducting phases have essentially the 
same characteristics. Thus it is hard to precisely con¬ 
clude the reason for the failure to detect such transition 
with our measure of entanglement. 


The discontinuities in the entanglement are directly 
identified with the first order quantum phase transitions, 
whereas the minimum points are identified with the sec¬ 
ond order quantum phase transitions. When crossing a 
first order transition, the ground state energy presents a 
discontinuity and consequently also its observables. In 
this way, the eigenvalues “A*,” of the single particle re¬ 
duced density matrix (Eq. ©), and the entanglement 
obtained from them, should present a discontinuity. The 
occurrence of the minimum points are due to the diver¬ 
gence of the correlation length when approaching the sec¬ 
ond order transitions. As described in the previous sec¬ 
tion, the eigenvalues “Aare given in momentum space 
by the Fourier t ran sform of the real space quadratures 
(Eq.(24l). In this way, if we are close to the 

transitions, such real space quadratures tend to become 
delocalised or spread out along the lattice, thus lead¬ 
ing to more localised eigenvalue distributions in momen¬ 
tum space, and consequently to smaller von Neumann 
entropies. It is worth remarking that such behavior is 
the opposite of the entanglement of modes, where the 
sites are maximally entangled at the second order tran- 
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FIG. 3: (Color online) Contour map for the entanglement 
of particles a E p ” as a function of the interaction terms V/t 
and U/t, in a system with L = 128 sites at half-filling. The 
entanglement behavior in the thermodynamic limit, L — ¥ oo, 
keeping fixed the filling n = N/L = 1, is qualitatively the 
same, with slight differences of the order of C?(10 -2 ) in its 
magnitude; see Appendix [X] for a detailed discussion. The 
(green) continuous line denotes the discontinuity at the en¬ 
tanglement function, while the (red) dashed line denotes the 
local minima. The white dots correspond to the points where 
we performed a detailed finite-size scaling analysis (see Table 
I in Appendix |A|) . 

sitions. As an example, see in Fig|4]the eigenvalue dis¬ 
tribution for a system with L = 128 sites when crossing 
the BOW-CDW quantum phase transition. 

We present now the entanglement behavior in some 
specific slices of the phase diagram with L = 128, in 
order to clarify the above discussion and results. More 
specifically, we show the entanglement behavior in the 
PS-SS-CDW, PS-SS-SDW, and PS-SDW-CDW transi¬ 
tions. Notice that our finite-size scaling analysis showed 
that in thermodynamic limit the entanglement behavior 
is qualitatively similar (see Appendix), with a scaling in¬ 
versely proportional to the lattice size, E p = aL -1 + b , 
where a and b are constants. 


A. PS-SS-CDW 

In Fig(5] we see the entanglement behavior across the 
PS-SS-CDW phases. We clearly see, for any fixed at¬ 
tractive on-site interaction (U/t < 0), a discontinuity in 
the entanglement followed by a local minimum point, as 
we increase the value of the inter-site interactions V/t. 
The discontinuity is related to the first order transition 
PS-SS, while the local minimum is related to the second 
order transition SS-CDW. We see, however, that the SS- 
CDW transition is not located exactly at V/t = 0, as 
expected from the phase diagram described in the liter¬ 
ature, but at a value close to this one. We believe that 



0 2 4 6 

k 

FIG. 4: (Color online) (top) Single particle quadratures 
at /2 «,V along the lattice sites, and (bottom) eigenvalue 
distribution “At,” for the single p arti cle reduced state in a 
fixed spin sector, as given in Eq.( |24[ ). We consider a fixed 
U/t = 4. The vertical axis is in log-scale, in order to make 
clearer the visualisation. As we approach the BOW-CDW 
quantum phase transition point, at V/t ~ 2.13, we see that 
the quadratures tend to delocalise along the lattice, whereas 
the eigenvalue distribution becomes more localised. 


this discrepancy is related to finite-size effects. 


B. PS-TS-SDW 

In Fig(6]we see the entanglement behavior across the 
PS-TS-SDW phases. We see again the two kinds of 
behavior for any fixed attractive inter-site interaction 
(V/t < 0): a first discontinuity, related to the first order 
transition PS-TS, followed by a local minimum point re¬ 
lated to the second order transition TS-SDW. Note that, 
for large values of the attractive inter-site interaction, 
V/t ~ —1.5, the discontinuity and minimum converge to 
the same point, and there is no TS phase anymore. 
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FIG. 5: (Color online) Entanglement behavior across the PS- 
SS-CDW phases. The entanglement, for any fixed atractive 
on-site interaction ( U/t ), is characterized by a discontinuity 
(PS-SS transition), followed by a local minimun (SS-CDW 
transition). 



FIG. 6: (Color online) Entanglement behavior across the PS- 
TS-SDW phases. The entanglement, for any fixed atractive 
inter-site interaction ( V/t ), is characterized by a discontinuity 
(PS-TS transition), followed by a local minimun (TS-SDW 
transition). For large V/t, the two transitions shrink at the 
same point, and there is no TS phase anymore. 

C. PS-SDW- (BOW)-CDW 

In FigjT] we see the entanglement behavior across the 
PS-SDW-BOW-CDW phases. We see that, as we in¬ 
crease the value of the inter-site interactions, for any 
fixed repulsive on-site interactions (U/t > 0), the en¬ 
tanglement identifies two transitions. Firstly we see a 
discontinuity, related to the first order transition PS- 
SDW, followed then by: (i) a discontinuity, when con¬ 
sidering large U/t, or (ii) a local minimum point, when 
considering small U/t. Such discontinuity is related to 
the first order SDW-CDW transition, while the minimum 
points are related to the second order BOW-CDW tran¬ 




FIG. 7: (Color online) Entanglement behavior across the 
PS-SDW-BOW-CDW phases. The bottom panel is a magni¬ 
fication of the top panel in the region 1 < V/t < 4.5. The en¬ 
tanglement, for any fixed repulsive on-site interaction (U/t), 
is characterized by a discontinuity (PS-SDW transition), fol¬ 
lowed by: (i) a discontinuity for large V/t (SDW-CDW transi¬ 
tion), or (ii) local minimun points for small V/t (BOW-CDW 
transition). 


sition (the SDW-BOW transition is not seen, as afore¬ 
mentioned). We see that the transitions to the CDW 
phase occur at U ss 2V. Performing a finite-size scaling 
analysis (see Appendix) we obtain that, for U/t = 4, the 
BOW-CDW transition is located at V/t = 2.11 ± 0.01, 
which is slightly lower than the literature results, namely 

v/ttt 2.i6 mmmm- 


V. CONCLUSION 

We studied the entanglement of indistinguishable par¬ 
ticles in the extended Hubbard model at half-filling, with 
focus on its behavior when crossing the quantum phase 
transitions. Our results showed that the entanglement 
either has discontinuities, or presents local minima, at 
the critical points. We identified the discontinuities as 
first order transitions, and the minima as second order 





































transitions. In this way, we concluded that the entangle¬ 
ment of particles can “detect” all transitions of the known 
diagram, except for the subtle transitions between the 
superconductor phases TS-SS, and the transition SDW- 
BOW. 

It is also interesting to compare our results with other 
entanglement measures, such as the entanglement of 
modes, which was widely studied in several models, as 
well as in the extended Hubbard model IHlTj . Gu et 
al. |38j firstly showed that the entanglement of modes, 
i.e., the entanglement of a single site with the rest of the 
lattice, could detect three main symmetry broken phases, 
more specifically, the CDW, SDW and PS. Other phases 
were not identified due to the fact that they are asso¬ 
ciated to off-diagonal long-range order. Further investi¬ 
gation were performed analysing the block-block entan¬ 
glement E3 EH, i.e., the entanglement of a block with 
l sites with the rest of the lattice (L — l sites), show¬ 
ing that this more general measure could then detect the 
transition to the superconducting phase, as well as the 
bond-order phase. The measure, however, could not de¬ 
tect the SS-TS transition, besides presenting some unde¬ 
sirable finite-size effects in the PS phase. On the other 
hand, the entanglement of particles studied in this work 
showed no undesirable finite-size effects in the PS phase, 
but could not detect the superconductor SS-TS transi¬ 
tion either. Regarding the BOW phase, from the above 
discussion we see that it would be worth to analyze more 
general measures for the entanglement of particles, which 
goes beyond single particle information. Some steps in 


this direction were made in [!T5] , where a notion of entan¬ 
glement of “subgroups” of indistinguishable particles was 
defined. 
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Appendix A: Finite-size scaling analysis 

In this appendix we perform a finite-size scaling anal¬ 
ysis in the system entanglement, in order to extract in¬ 
formation about the ground state of the model. We ob¬ 
tained that the entanglement behavior is qualitatively 
the same for lattices larger than L ss 100, with just small 
differences of the order of 0(1O -2 ) in its magnitude. In 
a general way, the entanglement scales with the inverse 
of the lattice size, E p = aL _1 + 6, where a and b are con¬ 
stants. See in Fig(8], for example, the entanglement scal¬ 
ing for the SDW-BOW-CDW phase transitions. In Tab. 
I we show the computed values for the scaling constants 
at different points in the phase diagram, as highlighted 
in Fig (5] 
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highlighted by the scaling at V/t = 4.147 (bottom curve), and 
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aL~ x + E P (L —> oo), at different points of the phase diagram, 
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cal points, and ” means that the entanglement is constant, 
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